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Abstract. We consider dissipative systems resulting from the Gaussian and alpha-stable noise 
perturbations of measure-preserving maps on the d dimensional torus. We study the dissipation 
time scale and its physical implications as the noise level e vanishes. 

We show that nonergodic maps give rise to an 0(l/e) dissipation time whereas ergodic toral 
automorphisms, including cat maps and their d-dimensional generalizations, have an 0(ln(l/e)) 
dissipation time with a constant related to the minimal, dimensionally averaged entropy among 
the automorphism's irreducible blocks. Our approach reduces the calculation of the dissipation 
time to a nonlinear, arithmetic optimization problem which is solved asymptotically by means of 
some fundamental theorems in theories of convexity, Diophantine approximation and arithmetic 
progression. We show that the same asymptotic can be reproduced by degenerate noises as well as 
mere coarse- graining. We also discuss the implication of the dissipation time in kinematic dynamo. 
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1. Introduction 

Irreversibility and approach to equilibrium are fundamental problems in statistical mechanics 
and dynamical systems and its complete solution is still elusive (see, e.g., JMI). There are possibly 
many routes to irreversibility. 

One view is that macroscopic systems are exceedingly difficult to isolate from their environments 
for a time comparable to their dynamical time scales. The noise as a result of interaction with 
environment may further trigger irreversibility, such as approach to equilibrium, in the systems. 
The initial uncertainty involved in preparing a physical system and the random perturbation due 
to measurements as well as Gibbs' coarse-graining procedure can all be viewed as certain noises. 
The point is that noises, intrinsic as a result of internal stochasticity or extrinsic as a result of 
random influence from surrounding environment, can induce effects that would be weak or absent 
without noise. 

In this paper we investigate one such effect, called dissipation, for discrete time, conservative 
dynamical systems under the influence of noise. In particular we study the time scale, called the 
dissipation time, on which the dissipation as measured in L*'— norm, 1 < p < oo, has an order one 
effect even as the magnitude of noise vanishes. Clearly the dissipation time depends on the ergodic 
properties of the noiseless dynamics as well as the noise level. 

The noisy dynamical system considered in this paper can be viewed as a discrete generalization 
of the dynamics of a passive scalar in a periodic, incompressible velocity field v 

(1) dx.%t) = v(x^(t))dt + ^/edw{t) 

V • v(x) = 0, 

where the standard Brownian motion w and the molecular diffusivity e represent the stochastic 
perturbations as a result of random molecular collisions (see, e.g., ^3]) |S])- The discrete-time 
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dynamical system will be defined on the d-dimensional torus T*^ = /lA. The velocity field v will 
be replaced by arbitrary Lebesgue-measure preserving map F defined on T'^ (periodicity condition). 

In order to study the dynamics generated by F it is useful to consider its Koopman operator Up 
defined by a composition Upf '■= f ° F, with / belonging to some Banach space of functions on 
T"^. We will be mainly concerned with the standard Banach spaces LP{T'^), 1 < p < oo, and their 
subspaces Lq{T'^) of functions with zero mean (/) = 0, where (/) denotes the average of / w.r.t. 
the Lebesgue measure. In case of L^(T"') one can consider Up as the Frobenius- Perron operator 
associated with F~^. 

In the time-discrete version we consider general a— stable noise operator Ge,a '■ -^oC^^) ^ -^o('^°')' 
with a G (0, 1], defined by means of the Fourier transform of corresponding a-stable noise kernel 

9e,a 

G,,,/(x)= / 5.,a(x-y)/(y)dy= ^e-W'"/(k)ek(x), 

where 



with ek(x) := e^'^''^'-^, k E Z'^. Here, just like in e > represents the level of the noise. Putting 
a = 1 one recovers standard heat kernel. 

The operator T^^a on Lq(T'^) generating the noise-perturbed dynamical system considered in this 
paper is thus given by 

(2) T,,^f:=G,,aUFf = ge,a*(.f°F)- 

Simple computations yield 

(3) \\T^,a\\ = \\Gs,aUF---Gs,aUF\\ < ||G£,q,||" = e~^". 

Here and throughout the paper || • || denotes the standard Lg-norm or the corresponding operator 
norm (any other norm will be equipped with suitable index). 

We define the dissipation time as the time on which the contraction © becomes of order one: 

(4) Udiss ■■= min{n G Z+ : ||r;;^|| < 1/e}. 

Hence the dissipation time is a function of e, a as well as the underlying dynamics. The choice 
of the threshold in the definition is a convenient one and for the purpose of the paper can 
be any positive number less than one (see Proposition P). The fact that for all e > 0, llT^all is 
monotonically decreasing ensures that ndiss is well defined. By contrast, when e = 0, the fine- 
grained Boltzmann-Gibbs entropy as well as L^-norm of the initial state remains constant in the 
course of evolution. In other words this is a "dissipation" effect and hence the term "dissipation 
time". On the dissipation time scale the system is, in a sense, "half way" through its irreversible 
route to the equilibrium state. The dissipation time provides a measure of the instability of the 
dynamics w.r.t. the stochastic perturbations which result in the "aging" of the system toward the 
final state. 

The main purpose of this paper is to investigate the asymptotics of the dissipation time as e 
tends to zero. Due to the non-normality of the operator T^^a, the dissipation time can not be 
determined from its spectral radius. Indeed, as we will see below, the operator T^^a corresponding 
to any ergodic toral automorphism F is quasi-nilpotent for any e > (thus the time scale estimated 
from the spectral radius is infinite) whereas the dissipation time is of the order In (1/e). 



To briefly describe the main results we pause to note the following asymptotic notation. Given 
two sequences , he indexed by the parameter e > we write 

Oe < 6e, if limsup 7^ < oo 
Op ~ 6p, if lim — ^ = 1 

Moreover we write ^ be, if both < be and be < hold simultaneously. 

The first result is that the dissipation time Udiss ~ 1/^ for nonergodic or, more general, non- 
weakly-mixing maps (cf. Theorem n and its corollaries. Section |2.3() . which is also the longest 
possible time scale for dissipation in view of © (see also Lemma ^ Section |2.2() . In other word, 
such systems are most stable w.r.t. stochastic perturbations. 

The main aim, however, is to investigate the cases in which the dissipation is much faster 
due to rapid mixing in F. We show (Theorem 2, Section 12. 3|) that, for a toral automorphism 
ndiss ~ log(l/£) if and only if the map F is ergodic (which in this case is also an Anosov 
difFeomorphism) . In particular our results hold for all classical cat maps (hyperbolic automorphisms 
of 2-torus) and their d-dimensional generalizations. In addition, we provide a general lower bound 
for the constant of the logarithmic asymptotics (Theorem 2). We further show that the lower bound 
is achieved for diagonalizable automorphisms, namely 

(5) Udiss « {2ah{F))-' ln(l/e) 

where h{F) denotes the minimal, dimensionally averaged entropy among F's irreducible blocks 
(Theorem 3, Section |2.3|) . Dimensionally averaged entropy for each irreducible sub-block of the 
toral automorphism is the Kolmogorov-Sinai (KS) entropy per dimension of an irreducible factor 
of the whole map. 

Our method involves solving asymptotically a quadratic arithmetic optimization (i.e. quadratic 
integer programming) problem by obtaining sharp upper and lower bounds using number theo- 
retical tools including multidimensional Diophantine approximation theorems (Schmidt's subspace 
theorem), Minkowski's theorem on linear forms and Van der Waerden's theorem on arithmetic 
progressions. This is done in Section |3J 

In Section IS] we show that the same result © holds when the noise is replaced by coarse-graining 
the initial and terminal states. This is reminiscent of the well-know results of statistical stability 
in the literature, namely, the Bernoulli systems are stable under the sufficiently small intrinsic 
random perturbation in the rough sense that the perturbed system is close to the direct product 
of the unperturbed one and some auxiliary viewer system (see jl3j . |21jl. In other words, for those 
systems, the process that results from small intrinsic random perturbation can be reproduced 
exactly by looking at the unperturbed system through a viewer that distorts randomly but slightly. 
In spite of the above the asymptotic Q indicates the perturbed system is irreversibly far from the 
unperturbed one even on the relatively short dissipation time scale. From this perspective, such a 
system is statistically unstable. 

In Section |1] we consider a class of highly degenerate noises and show that the same conclusions 
about the dissipation hold if the degenerate noises satisfy an additional generic condition. 

In Section El we consider the relation between the dissipation time and some characteristic time 
scales relevant to kinematic dynamo. We show that fast dissipation generally inhibits dynamo 
action. When there is no fast dynamo action the noisy push-forward map dissipates the magnetic 
field energy on the dissipation time scale. However, the magnetic field energy can still grow to 
relatively large magnitude as inverse power-law of the small noise with the exponent proportional to 
the ratio of the logarithmic spectral radius of the toral automorphism to the minimal, dimensionally 
averaged entropy among the automorphism's irreducible blocks (cf. (|65])1. 

The notion of dissipation time has a natural bearing on the problems of quantum chaos with noise. 
The family of symplectic toral automorphisms constitute important examples of quantizable chaotic 
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dynamics on compact manifolds for which various quantization procedures have been intensively 
studied (see, for example, i22j and |19j). We will address the issue of decoherence time for quantized 
symplectic toral automorphisms with noise in a forthcoming paper. 

The organization of the rest of the paper is as follows. In Section |21 we develop the general 
theory of dissipation time and its relation to the Boltzmann-Gibbs entropy. We also formulate the 
dissipation time calculation for total automorphisms as an arithmetic minimization problem and 
state the main results. In Appendix we generalize the dissipation time asymptotic result to the 
affine transformations. The proofs of some elementary facts are presented in Appendix^ for the 
sake of completeness. 



2. Dissipation time 

In its general form the dissipation time rifiigsip) can be defined in terms of the norm || • ||p^o on 
the space Lq(T'^) w.r.t. a threshold rj G (0, 1) 

(6) nrfj5s(p,r/) := min{n e Z+ : ||r"„||p,0 < 1 < p < oo. 

First we show that the value of the threshold in @ does not affect the order of divergence of 
f^diss{p-,v)i ^ ^ tends to zero. 

Proposition 1. For any < J?,?] < 1, ndiss{p,'n) ~ "-diss (p, ??) • 

Proof. Assume < fj < t] < 1. Obviously ndiss{p,fj) > ndiss{p,ri). On the other hand let k be 
a positive integer such that 7]^ < fj. Then 

(7) rj^^-^P'^^llp.o <r]^ WT^l'-^^^Wpfl <r,^<fl. 
Hence kndissip,v) > "^dissiPi v) : which implies n^iss ip,v) 

{P,fj). ■ 

Following the argument of |2SI one can use the Riesz convexity theorem to establish also the 
asymptotic equivalence of the Udissip), for all 1 < p < cxd. 

Proposition 2. i) For any 1 < q,p < oo, Udissiq) ~ n-dissip)- 

a) For any I < p < oo, Udtssip) < ridtssi'^) and Udissip) < ?^diss(oo). 

The details of the proof can be found in Appendix IbI 

Our particular choice of the exponent p = 2 and threshold rj = in Q) is computationally 
convenient and will be used throughout the paper. We will use the convention that Udissip) = 
ridissip, e~^). 

We say that operator T^^a or associated with it measure preserving map F has a simple (slow) 
dissipation time when Udiss ~ 1/e and that it has a logarithmic (fast) dissipation time when 
Udiss ~ ln(l/e). 

In the particular case of fast dissipation, with a logarithmic dissipation time, in order to estimate 
precisely the rate of dissipation, one needs to determine the value of the dissipation rate constant 
Rdissi defined as 

(8) Rdiss ■= lim "'^'^^ 



oln(l/e) 

Similarly in case of simple dissipation time the dissipation rate constant can be defined as 
(9) Rdiss ■■= li'oiendiss- 

s—>0 



2.1. Dissipation time and Boltzmann-Gibbs entropy. In this section we briefly discuss the 
connection between dissipation time and Boltzmann-Gibbs entropy. 

First we note that on the scale of n^j^s the Boltzmann-Gibbs entropy approaches the maximal 
equilibrium value (i.e. 0) as can be seen from the following simple estimate |16| . Let us first restrict 
considerations to bounded initial states, i.e., / > 0, / G and ||/||i = 1. Let 



7]{u) 



—ulnu, u> 
0, u = 



and let Z?„ = {x G T*^ : 1 < T^^f}. On one hand, we have 



r?(r" /(x))dx 



< 



Dn 



dr]{u) 
du 



du 



dx 



< sup (l+lnu) [ |r,"^/(x) - l|dx 

1<«<||T" /lloo JD„ 



J\\ 

< (i + in||r,viloo)||r,V-i||i 
(10) < (i + in||/|U)||r,v-i||i- 

On the other hand, we have 

> / r?(Tj:„/(x))dx > / r,(T-,/(x))dx. 

In view of the inclusion relation: L°°(T'^) C ^^(T'^) c ^^(T'^), we then obtain that for n ^ n^j, 

e[0 



sup 

/>0,||/|U<c 



r^{T^^J{^))d^ 



0, Vc > 0. 



For unbounded initial states, we note that, by Young's inequality, 

l|-^e,a/l|c!0 ^ ||?e,a/||oo ^ ||5e,a II oo || / 1| 1 — Hffe, alloc 

from which we have, instead of (|1U() . the following estimate 



r?(T" /(x))dx 



< (1 + In lift 



7;v-iiii. 



where 

In lls'e.alloo ~ ln(l/e). 
Therefore for sufficiently fast diverging n ^ ?^djss(l) such that 



(11) 

one obtains 



ln(l/e)||r,':J/-l)||i,o 

r?(r" /(x))(ix 







sup 

/>0,||/||i=l 



fd 



0. 



The condition (lllj) typically results in a slightly longer time scale than Udissi^)- 

On the other hand, we can bound the Li distance between the probability density function / 
and the Lebesgue measure by their relative entropy via Csiszar's inequality 



Yd 



|/(x) - 5(x)|<ix < J2 /(x) In (/(x)/<7(x))dx 



with ^(x) = 1. We see immediately that the decay rate of 



sup 

/>0,||/||i=l 



??(r- /(x))dx 



provides an estimate for Udissi^) and, consequently, for ndissip)^? £ (li oo). 

2.2. Calculating the dissipation time. For greater generality and transparency of arguments we 
consider, in this section, a slightly more general family of operators T^^q, defined, as previously, by 
the first equality in but with arbitrary unitary or isometric (not necessary Koopman) operator 
U (and hence in these cases we drop the subscript F). 

Lemma 1. For any isometric operator U , the dissipation time ofT^^a satisfies following constraints 
(12) \\R{l;T,,^)\\<naiss<l/e, 
where i?(l;T£.Q,) denotes the resolvent ofT^^a o,t 1. 
Proof. In view of Q and ©, for n = ridiss one has 

g-l < ||2^(«diss-l) II < g-e("disa-l)^ 

which clearly implies the second estimate of (fT^ . In order to prove the other inequality we proceed 
as follows. 

no — 1 oo 



no-1 



/ > £,0 



n=0 



n=0 



n=0 



Erpn I rpriQ \ ^ rpn 

n=0 n=0 

<no+||r"«||||i?(l;r,,,)||. 



Hence taking in the above inequality no = ridiss one gets 

||i?(l;r,,,)||(l - e'') < p(l;r,,„)||(l - ||r,^--||) < n,,,„ 

which gives the first estimate of ()12() . ■ 

The above lemma provides an absolute upper bound for dissipation time. Taking F = I one 
easily finds that this bound is best possible in general. The lower bound is useful in the case when 
one can estimate from below the norm of the resolvent (see proof of Theorem 

Theorem 1. If U acting on Lq(T'^) possesses nonempty pure point spectrum and at least one of 
its eigenf unctions belongs to H'^"{T'^), then T^^^ has simple dissipation time. 

Proof. In view of Lemma ^ it is enough to find a lower bound for the norm of the resolvent 
i?(l;T£^Q,). Let h G H^"' be one of the eigenfunctions of U. Since U is isometric we have 

Uh = e'^h. 



We first assume that 

\\m;Te,a)\\ 



0. Since 1 ^ a{T£^a), I — T^^a is a homeomorphism and hence 



sup 



il - T,, 



sup 



> 



^,rg||(/-r,,„)/|| - 11(1 -r,,,) 



Now expressing h in the Fourier series we get 



(i-r,,„)/if = ^ |/i(k)p 



1 



-e|k|- 



< (^IMk)llkl 



2a 



2 

^2o 



Hence 



> 



e\\h\ 



C 

e 



Thus in view of (|12|) and above calculations 

l/e<\\R{l-Te,a)\\<ndiss<l/e, 

which ends the proof in the case (/> = 0. 

If (/) 7^ 0, we put {/ = e~^^U, which implies JJh = h. 

The proof is completed by applying the above reasoning to operator q, = G^^aU and observing 
that the dissipation times for T^^a and Tg^Q, are identical. ■ 

When [/ is a Koopman operator associated with a map F, then the property that U considered 
on Lq(T'^) possesses nonempty pure point spectrum is equivalent to the fact that F is not weakly 
mixing (see (SI). Thus we have 

Corollary 1. If F is not weakly mixing and its Koopman operator possesses H^°' eigenfunction in 
Lq{T'^), then T^^^ has simple dissipation time. 

Another immediate consequence is 

Corollary 2. If F is not ergodic and its nontrivial invariant measure possesses H^" density func- 
tion, then Ts^a has simple dissipation time. 

A typical example of ergodic but not weakly mixing transformations for which the above corollary 
applies is the family of 'irrational' shifts on T'^ i.e. maps Fx = x + c on T"^, where c = (ci, .., c^) is 
a constant vector such that the numbers l,ci, ..,Crf are linearly independent over rationals. More 
general and less trivial examples of ergodic maps giving rise to a simple dissipation time will be 
discussed in Appendix A (cf. Remark 

In general the problem of computing the dissipation time is rather complicated. In some cases 
it can be reformulated as an asymptotic optimization problem. To see it, one can represent the 
action of a given unitary operator U in the Fourier basis 



where for each k 
(13) 

Next we introduce the notation 
Z^„(ko,k„) = 

5„(k„) = 

Then for any / G Lq(T'^) we have 



1. 



7/ko,ki.-^ik„_i,k„e 

07^ki,...,k„_iGZ<* 

{koGZ^{0}:^/„(ko,k„)/0}. 



J2 /(ko)Tj:„eko 

O^koGZ'* 



(14) 



Yl Yl ^n(ko,k„)ek„ 

O^koGZ'* OT^knGZ-i 
2 



Y /(ko)^n(ko,k„ 

OT^koGZ'* 



E 

O^knGZ'* 



0^k„GZ'' 

The following general upper bound for ||TJ^q/II holds. 
Lemma 2. For any f G Lq(T'^), 

(15) WT^^jf < E E l/>o)l' E K(ko,k„ 



E /(ko)Z^n(ko,k„ 

koG5„(k„) 



0^k„GZ'* koG5n(k„) 



koG5„(k„) 



For the proof we refer the reader to Appendix IbI 
When «k,k' is a Kronecker's delta function 

(16) Uk,k' = ^Ak.k', 

where A : Z'^ i-^ Z'^ is a linear surjective map, the upper bound (|15|) can be used to obtain an 
identity for ||r"Q,||- First observe that 

and hence (|15|) becomes 



\KJf < E |/(ko)|V2^^-'^'''°'"<||/f max e-2^S-l^''^l^^ 
On the other hand for any nonzero k G Z"^, one can take in 1)14^ / = ek and get 

and therefore 

(17) ||r,"J| = max e-^S"-l^''^l'"=e-^'^^°o^''«'*EL.|A'k|-_ 

Let us now determine the class of maps F such that the corresponding Koopman operator Up 
satisfies (|16|). The relation (|16j) implies 

^7Fek = e^k = e2-^<^''''^>. 

On the other hand 

t/^ek(x)=ek(Fx) = e2-<'^'^'^>. 

Thus 

(k, Fx) = {Ak, x) mod 1, Vx G M"^, k G Z*^, 

that is, A is linear and A^ equals the lifting of F from T'^ onto M.'^. Moreover, the matrix A has 
integer entries and determinant equal to ±1, i.e., A (and F) is a toral automorphism. Hence, for 
toral automorphisms, the calculation of the dissipation time reduces to the following nonlinear, 
asymptotic (large n) arithmetic minimization problem 

n 

(18) min yi^^kp". 

We will show in Section |21 that for any ergodic toral automorphism this minimum value grows 
geometrically in n with the base related to the dimensionally-averaged KS-entropy of the total 
automorphism. 

2.3. Dissipation time of toral automorphisms. It is well known (see [l]) that (the lifting map 
corresponding to) any toral homeomorphism H : T'^ T'^ can be decomposed into three parts 
H = L + P + c, where L, the linear part, is an element of SL[d,'L) - the set of all matrices with 
integer entries and determinant equal to ±1, P is periodic i.e. P(x + v) = P(x) for any v G Z*^, 
and c is a constant shift vector. 

Every algebraic and measurable automorphism of the torus is continuous. Each continuous toral 
automorphism is a homeomorphism with zero periodic and constant parts and hence can be iden- 
tified with an element of SL{d,'L). And vice versa, each element of SL{d,'L) uniquely determines 
a measurable, algebraic toral automorphism. Thus from now on the term toral automorphism will 
simply be reserved for elements of SL{d, Z). We recall here that all Anosov diffeomorphisms on T*^ 
are topologically conjugate to the toral automorphisms ([Hj, [18]). 

Below we summarize some ergodic properties of toral automorphisms (cf. ^) p. 160, 12 and 
0). 



Proposition 3. Let F be a toral automorphism. The following statements are equivalent 

a) no root of unity is an eigenvalue of F. 

b) F is ergodic. 

c) F is mixing. 

d) F is a K-system. 

e) F is a Bernoulli system. 

In the sequel we will use the following result (cf. j26|). 
Proposition 4. The entropy h{F) of any toral endomorphism F is computed by the formula 



where Xj denote the eigenvalues of A. 

^From the formula (|19|) one immediately sees that a toral automorphism has zero entropy if and 
only if all its eigenvalues are of modulus 1. In fact much stronger result holds. 

Proposition 5. A toral automorphism has zero entropy if and only if all its eigenvalues are roots 
of unity. In particular all ergodic toral automorphisms have positive entropy. 

Given any toral automorphism F we denote by P its characteristic polynomial and by {Pi, P^} 
the complete set of its distinct irreducible (over Q) factors. Let dj denote the degree of polynomial 
Pj and hj the KS-entropy of any toral automorphism with the characteristic polynomial Pj. For 
each Pj we define its dimensionally averaged KS-entropy as 



Now we state two main theorems of the present paper. 

Theorem 2. Let F be any toral automorphism, Up the Koopman operator associated with F , G^^a 
a-stable noise operator and T^^q = Gs^aUp- Then 

i) T^^a has simple dissipation time if and only if F is not ergodic. 

a) T^^a has logarithmic dissipation time if and only if F is ergodic. 

Hi) IfTe^a has logarithmic dissipation time then the dissipation rate constant satisfies the follow- 
ing constraint 



Part i) of the above theorem follows immediately from Theorem ^ For details of a simple proof 
we refer to appendix IbI 

The natural question arises, whether the lower bound for the dissipation rate constant given in 
the above theorem is best possible. The next theorem and its corollary provides a strong argument 
in favor of this conjecture. 

Theorem 3. // F is ergodic and diagonalizable then 



(19) 




(20) 




For the whole matrix F we define its minimal dimensionally averaged entropy (denoted h{F)) as 

h(F) = min hj 
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That is, the dissipation rate constant ofT^^a is given by 

1 



R, 



diss 



2ah{F) 

The proof of parts ii) and iii) of Theorem [^J and of Theorem |21 constitute the most important 
part of this work and will be presented in Section 13.31 after necessary tools are developed. 

We end this section with the the results for two and three dimensional tori. Ergodicity of two 
dimensional toral automorphisms is equivalent to hyperbolicity. Two dimensional hyperbolic toral 
automorphisms are often referred to as the cat maps. 

Using Corollary |1] and applying Theorem |31 to two and three dimensions one gets the following 

Corollary 3. Let F he any ergodic, two or three dimensional toral automorphism. Then 

ndiss ~ — TTT 
2ah{F) 

3. Asymptotic arithmetic minimization problem 

In this section we find the asymptotics, as n goes to infinity, of the following quadratic arithmetic 
minimization problem 



(21) min yiA'kP", 

where A E SL{d,'Z). When A is not ergodic the asymptotics of ()21|) is of the order 0{n). For the 
rest of the paper we will consider only the ergodic case. For d = 2 the problem 1)21(1 can be solved 
easily as follows. Consider first the case that A is symmetric and a = 1. ^From det{A) = 1 we see 
that eigenvalues are A,A~^ with |A| > 1. We have 

2n+l n 



mm 



y i^'kp = min y l^'kl 



2 



(=1 l=—n 



min |k|2 + y |A|2'|ki|2 + |A|-2'|k2|2 + V \\\-^'\ki\^ + \\f\k2\ 



2 



n n 

21 



= min y |A|2'|k|2 = y |A| 

l=—n l=—n 

Hence there exist constants Ci and C2 such that 



Cie'^^^)" < min V M^kP" < C2e 



h{A)n 



where h(A) denotes the KS-entropy of A. The estimates for the general case of non-symmetric A 
and Q 7^ 1 are similar. 

In higher dimensions, the solution to H21() is much more involved because of the presence of 
different eigenvalues with absolute values bigger than one. We have the following general estimate 

Theorem 4. Let A E SL{d,'L) he ergodic. There exist constants Ci and C2 such that for any 
Q < 5 <1 and sufficiently large n 



(22) ^^g(l-5)2ah(A)n < ^-^^ |^Zj^|2a < (j^^^2ah{A)n 



where as before h{A) denotes minimal dimensionally averaged entropy of A and h{A) denotes a 
constant satisfying < h{A) < h{A), with equality achieved for all diagonalizable matrices A. 
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The question whether the equality h{A) = h{A) holds for all ergodic matrices remains open. 
The proof of the theorem relies on nontrivial use of three number-theoretical results stated below. 

I. Minkowski's Theorem on linear forms 

Let Li,...,Ld be linearly independent linear forms on which are real or occur in conjugate 
complex pairs. Suppose oi, 02, are real positive numbers satisfying 0102- --0^ = 1 and ai = Oj, 
whenever Li = Lj. Then there exists a nonzero integer vector k G such that for every j = 1, ...,d, 

(23) \Ljk\ < Dttj, 

where D = | det[Li, L^]] V-^. 

Minkowski's Theorem on linear forms will be used to obtain a sharp upper bound on the as- 
ymptotic solution of the arithmetic minimization problem. The proof of the above theorem and its 
generalization to arbitrary lattices can be found in j2Uj (Chap. VI). 

II. Schmidt's Subspace Theorem 

Let Li,...,L(i be linearly independent linear forms on with real or complex algebraic coeffi- 
cients. Given 6 > 0, there are finitely many proper rational subspaces ofW^ such that every nonzero 
integer vector k with 

d 

(24) JJ |Ljk| < |kr^ 

i=i 

lies in one of these subspaces. 

Schmidt's Subspace Theorem will be used in conjunction with Van der Waerden's Theorem on 
arithmetic progressions (see below) to obtain a sharp lower bound for the asymptotic solution of 
the arithmetic minimization problem. The proof of Schmidt's Subspace Theorem can be found in 
[21j (Theorem IF, p. 153). 

Definition 1. For a given set of linear forms and for fixed (5 > 0, the smallest collection of proper 
rational subspaces of which contain all nonzero integer vectors satisfying i24\ ), is called the 
exceptional set and denoted by Eg. 

A main difficulty to be resolved in using Schmidt's Subspace Theorem is to show that the 
minimizer of either the original problem (|21() or an equivalent problem does not lie in the respective 
exceptional set which is in general unknown. We will pursue the latter route by using Van der 
Waerden's Theorem on arithmetic progressions to show that one can always construct an equivalent 
minimization problem whose minimizer is guaranteed to lie outside the corresponding exceptional 
set. To this end we note that Schmidt's Subspace Theorem is true when the standard lattice 
Z*^ is replaced by any other rational lattice, that is any lattice of the form A = Q(Z'^) where 
Q £ GL(d, Q). Schmidt's subspace theorem can be generalized to this situation by considering the 
set of new forms Lj = LjQ. The fact that Q G GL{d, Q) implies immediately that Lj are still 
linearly independent forms on with real or complex algebraic coefficients. 

III. Van der Waerden's Theorem on arithmetic progressions 

Let k and d be two arbitrary natural numbers. Then there exists a natural number n^[k.,d) such 
that, if an arbitrary segment of length n > n^, of the sequence of natural numbers is divided in any 
manner into k (finite) subsequences, then an arithmetic progression of length d appears in at least 
one of these subsequences. 

The original proof was published in j^; Lukomskaya's simplification can be found in |17j . 

Before presenting the proof of our main results we state a number of technical facts concerning 
the structure of toral automorphisms. 

3.1. Algebraic structure of toral automorphisms. In this section we denote by GL{d,Q) 
the group of nonsingular d x d matrices with rational entries or the group of linear operators on 
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Euclidean space W^, which are represented in standard basis by such matrices. We generally use 
the same symbol to denote both operator and its matrix. 

In the sequel a vector x S will be called an integer (or integral) vector if all its components 
are integers, and similarly a rational, an algebraic vector if all its components are rational or 
respectively algebraic numbers. The term rational subspace ofW^ will then refer to a linear subspace 
of spanned by rational vectors (cf. p. 113). 

Definition 2. A € GL{d,i 
ducihle in Ofxl. 



is called irreducible (over Q) if its characteristic polynomial is irre- 



Lemma 3. The following statements about a matrix A G GL{d,Q) are equivalent. 

a) A is irreducible. 

b) A does not possess any proper rational A-invariant subspaces ofW^. 

c) No rational proper subspace ofW^ is contained in any proper A-invariant subspace o/M'^. 

d) For any nonzero q G Q'^ and any arithmetic progression of integer numbers ni, .■.,nd, the set 
{^"■iq, A"2q, ^"''q} forms a basis o/M'^. 

e) a'' is irreducible. 

f) No nonzero q E Q"' is orthogonal to any proper A-invariant subspace o/M°'. 

g) No proper A-invariant subspace ofW^ is contained in any proper rational subspace o/M'^. 

Definition 3. We say that operator A £ GL{d, Q) is completely decomposable over Q if there exists 
a rational basis of M*^ in which A admits the following block diagonal form 



(25) 



where for each j = 1, r < d, Aj G GL(dj,Q) is irreducible and Y2'j=i ~ 

In general, any matrix A G GL[d, Q) admits a rational block diagonal representation [Aj]j=i^...^r. 
The smallest rational blocks to which A can be decomposed are called elementary divisor blocks. 
The characteristic polynomial corresponding to any elementary divisor block is of the form p*", 
where p is an irreducible (over Q) polynomial (see, e.g., [7]). Although elementary divisor blocks 
cannot be decomposed over Q into smaller invariant blocks, some elementary divisor blocks may 
not be irreducible. This happens iff m > 1 iff A is not completely decomposable over Q. One has 
the following elementary fact (see Appendix^ for a proof). 

Proposition 6. A £ GL{d, Q) is completely decomposable over Q iff A is diagonalizable. 

However, even if A G GL(d, Q) is not completely decomposable, each elementary divisor block 
of A can be uniquely represented (in a rational basis) in the following block upper triangular form 

\B C 
D ' 





. 


. 







. 





. 





(26) 



where B is the unique rational irreducible sub-block associated with A-invariant rational subspace 
of that elementary divisor and C, D denote some rational matrices. 

Proposition 7. All the eigenvalues of an irreducible matrix A G GL{d,(Q) are distinct (complex) 
algebraic numbers. In particular all irreducible matrices are diagonalizable. 

The proofs of the above propositions can be found in Appendix B. 

Finally we note that since the leading coefficient and constant term of a characteristic polynomial 
of any toral automorphism are equal to 1, the only possible rational eigenvalues of such map are 
±1 or zizi. The latter fact implies that ergodic toral automorphisms do not possesses rational 
eigenvalues. Thus we have the following 
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Corollary 4. Let F be an ergodic, two or three dimensional toral automorphism. Then F is 
irreducible (and hence diagonalizable). 



3.2. Proof of Theorem |1J This section is entirely devoted to the proof of Theorem |11 

Let [^j]j=i,...,r be a rational block-diagonal decomposition of A into elementary divisor blocks. 
Since A G SL{d,'£), there exist a transition matrix Q G SL(d,Q) such that for every / G Z, 

(27) A^ = Q-H[A^]yQ 

and moreover each elementary divisor block [A]j is represented in its block upper triangular form 
pel). 

The matrix Q defines a new lattice A = Q{'L'^) and acts bijectively between this lattice and the 
standard lattice TJ^. Hence 

n n n 

Moreover 

||Q||-2°|([yl,])'q|2- < |Q-i([yl,])'q|2" < f"] ([A,]yq|2", V/,j,a. 

Now we decompose A into the direct sum of lower dimensional sublattices Aj corresponding to 
invariant blocks \Aj\. So that 



(28) min V|([Al)'q|2° = min min V|(A)'q|2° 



Thus, without loss of generality, we may specialize to the case that A is already indecomposable 
over Q i.e. A does not possesses any proper elementary divisor blocks. To simplify the notation 
we will work with the standard lattice A = Z"^. According to the remarks following the statements 
of Minkowski's and Schmidt's Theorems the proof can be easily adapted for any rational lattice 
K = Q{Z'^). 

Since the technique of the proof differs depending on diagonalizability of A we consider two cases: 

3.2.1. Diagonalizable case. Here we concentrate on the case when A is diagonalizable and hence 
due to its in-decomposability irreducible (cf. Proposition inj . 

We denote by Aj (j = 1, d) the eigenvalues of A. Following Proposition d we note that Xj are 
distinct (possibly complex) algebraic numbers and hence there exists a basis (of C^) {vj}j=i^...^d 
composed of normalized algebraic eigenvectors corresponding to eigenvalues Xj. 

We denote by [Pj]^^^ the projections on [vj], and by [Lj\ the corresponding linear forms. It is 
easy to check that [Lj\ are given, in the Riesz identification, by the eigenvectors of the matrix 
which are co-orthogonal to [vj], i.e., (uj,Vj) = for i ^ j. [uj] and [vj] are real or occur in 
complex conjugate pairs. We have 

d d d 

X = ^ P,x = = Y.{^, u,-)v„ Vx G R'^. 

j=i j=i j=i 

The equivalence between any two norms in a finite dimensional vector space, implies the existence 
of absolute constants Ci , C2 such that 

d d 

Cij;|P,-x|2<|x|2<C2 5;|P,-xp. 
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Using the above inequalities, the monotonicity of a map x i-^ x" and an obvious inequahty (a+6)" < 
a" + 6", which holds for all positive o, b and a E (0, 1] one obtains the following estimates 

n n I d \ n / d ^ 

1=1 1=1 \ j=i j [=\ \j=i 

n d d / n \ 

1=1 j=i j=i \i=i / 

and on the other hand 

1=1 \i=i / \i=i j=i 

n d \ d 



ci I Y.Y.\H''\m = E Ei^/' i^M 

=1 i=i / Vj=i \i=i / 



2 



Now we introduce some notation 

(29) := max{l,|Aj|}, 

. d ^'1' 

(30) Vo := (HA,- 
One can easily observe that there exists a constant C such that 

n 

CAf" < J];|Ajf°' < nAf". 

In the sequel we do not distinguish between particular values of constants appearing in computa- 
tions. The symbols Ci, C2, .. are used to denote any generic constants independent of n. 
The normalization condition |vj| = 1 implies the following relation 

(31) |Pjx| = |Ljx|. 
Combining the above estimates one gets the following general bounds 

(d \ n d 

^Af |Ljk|2 < ^|A'k|2" < C72n^Af"|Ljk|2". 
j=i / 1=1 j=i 

Therefore in order to estimate (PT|l it suffices, essentially, to estimate 

d 

(33) min V A?" "iL.kP". 

We denote by z„ the sequence of minimizers i.e. nonzero integral vectors solving H33|l . 

Upper bound. For the upper bound we assign to the set of linear forms Lj the set A composed 
of all real vectors a = (ai,...,Orf) satisfying the conditions Oj > 0, for j = and Oj = 

whenever Lj = Lj and 

d 

(34) n«i = l- 
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/,From Minkowski's theorem on linear forms, we know that for any a € A, there exists nonzero 
integral vector ka satisfying |Ljka| < Daj, j = 1, d, where D = \ det[Li, L^W^^^. 
Thus 

d d 

(35) Af "iLjkap" ^ ^ • 

The minimizing property of z„ implies that for any a G 

d d 

(36) Yl "|LjZ„|2° < Y Af "|L,ka|2". 

i=i i=i 

Thus combining (|35jl and (|36|1 . and applying the Lagrange multipliers minimization with the con- 
straint (|34j) (and using the fact that Aj = Xj whenever Li = Lj), we get 

d d ( V'^ 

(37) Y Af "l^jZnl^" < D min ^ Af "af = '^^11 = '^^^^ 



i=i i=i \i=i 



Thus the following upper bound holds 

n 

(38) min V |A'k|2° < CanA^^^". 

Lower bound. Let m denote an arbitrary natural number. Using the fact that A acts bijectively 
on 'L'^ we can restate the minimization problem (|33|) in the following form 

d d 

(39) min V A?""|L,kP° = min V A?""|L,.4-™A'"k| 



2a 



mi,| 2o 



1 2a 



(40) = min y A?"" I A, 1-2°™ I L. A 
That is 

d d 

(41) ^ Af "|LjZ„|2" = ^ Af "|Aj-|-2""|Lj-A" 

i=i i=i 

We choose arbitrary 5 > and consider the exceptional set Es (see Definition^ associated with 
the system of linear forms [Lj]. Since [Lj] correspond to the eigen-pairs [Aj,Uj] of A'^ they are 
linearly independent linear forms with (real or complex) algebraic coefficients. Thus the subspace 
theorem asserts that Es is a finite collection of proper rational subspaces of W^. We denote by ks 
the number of subspaces forming Es- 

Now we want to show that for all sufficiently large n there exist an integer m < n such that 
A'^Zn does not lie in any element of Es- To this end we assume to the contrary that all A"^Zn lie 
in the subspaces forming Eg and we divide the sequence of natural numbers 1, ...,n into kg classes 
in such a way that two numbers mi and 1712 are in the same class if A^^Zn and A'^^Zn lie in the 
same element of Eg. Now let ni,{ks, d) be the number given in the van der Waerden theorem and 
let n > n*. Then there exists an arithmetic progression mi,...,md in one of these subsequences. 
By LemmaEld) the set of vectors {A^^Zn, A^'^Zn, ...,^"*<*z„} forms a basis of the whole space M*^, 
which contradicts the fact that they lie in one fixed rational proper subspace. Hence for any 5 > 
and n > riif there exists < n such that A'^'Zn does not lie in any element of Eg. 

Now, introducing the notation 

(42) in = A'^'Zn 
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one concludes from H41() that for any S > and all n > the following equality and estimate hold 

d d 

(43) ^Af"|LjZ„|2'^ = ^Af"|Aj|-2"'^*|LjZ„|2" 

d 

(44) ni^j-^"i ^ 

Inequality may be rewritten as 

d 

(45) n 1^^-^ 
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with some / : — > such that /(r) < r,\/r > 0. 

Using (|^^ and (|H7|l we obtain the existence of a constant A > 1 such that 

d 

(46) /(|Z„|) < |Z„| = |^™*Z„| < A:^;Jz„| < X^,, Yl >^l\Lj^n\ < dD(\maAgeoT < A"- 

Note that Hj Xj = 1. So, by (g^) the quantities 5^,^ = (|Aj|-"*/(|z„|)'5/'^|LjZ„|)^" , j = l,...,ci 
satisfy the constraint 

d 

(47) J]5j.„ = l, Vn>n,. 

i=i 

Thus applying (|46|) and the Lagrange multipliers minimization with the constraint (|47jl one gets 

d d 

Xj \Xj\ \LjZn\ =/(|Zn|j ' 2^Xj Bj^n>X ' Xg^^ =. Xg^^^ '. 

This and equality (|43|) yields the following lower bound for (|33|) 

(48) min V Af "iL.kP" > ^""(i-^). 

3.2.2. Non-diagonalizable case. We move on to the general case where A is not irreducible (but, 
as assumed at the beginning of the proof, indecomposable over Q). We denote by B the invariant 
irreducible sub-block of A given by its block upper triangular decomposition ()26p and by S the 
rational invariant subspace associated with this block. We note that B as an irreducible matrix is 
diagonalizable. 

Upper bound. Note that 

n n 

(49) min V M'kP° < min V l^'kP". 

The corresponding upper bound H38|) for i3 is then also an upper bound for the whole matrix A. 
We note that geometric average of Xj over S is equal to the geometric average of all Xj associated 
with matrix A (i.e. over the whole space M'^). 

Lower bound. According to our assumption A is indecomposable and thus the characteristic 
polynomial of A is of the form p"^ for some irreducible p. All Jordan blocks of A have the same 
size m and different Jordan blocks correspond to distinct eigenvalues. Denote by b the number of 
the Jordan blocks in A and by Xj, where j = 1, all these distinct eigenvalues . Since each Xj 
has algebraic multiplicity m, we get d = mb. Let {"Vj h}j=i^ h-h=o,...,m-i be a basis (of C^) in which 
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A admits the Jordan canonical form. As usually Lj^h will denote the corresponding linear forms. 
Each Vj can be regarded as a generalized eigenvector of A associated with an eigenvalue Xj. We 
assume that these generalized eigenvectors are ordered according to their degree i.e. Vj^h satisfies 
the equation {A — XjlY^^Wj^h = 0. Reordering the eigenvalues, if necessary, we can also assume 
that Ai has the largest modulus among all eigenvalues of A and hence Ai = |Ai|. Let z„ be the 
sequence of minimizers solving (|21() . We first note that for each n there exists < /i < m — 1 such 
that Li /jZ„ 7^ 0. Indeed, otherwise for all /i = 0, ...,m — 1, Li /jZ„ = and consequently for any n 
and h Li ^A^Zn = 0. The latter implies that the set of consecutive iterations {zn, A^Zn, A'^Zn, ...} 
spans a proper rational A-invariant subspace of M'^ which does not have any intersection with the 
subspace spanned by the generalized eigenvectors of A associated with eigenvalue Ai. This clearly 
contradicts the irreducibility of p. Now, for given n we denote by h{n) the biggest index h for 
which the condition LijiZ„ ^ holds. 
We have the following estimate 



(50) 



(51) 



\ 2an I T „ 



< 




1-h 



E 



n 



L 



'j,h 



where the last inequality follows from previously established upper bound. 

/,From the Diophantine approximation and the assumption that |Li /i(„)Z^ 
/? > such that (see [21 p. 164) 

1 



7^ 0, there exists 



(52) |Li^/j(„)Z„ 

Thus combining (|5fl]) with ([52]) one gets 

A? 



> 



1/3' 



2«ni_ \—2a!3 ^ \2cm\T 



2a ^ 1^ ^\ 2an 



After rearrangements one obtains the following lower bound estimate for (|'21|) 



(53) 



where 



(54) 



-A 

n 



2an 



< C\z, 



\2a 



n 



< 2^ l^'zn 
1=1 



2a 



mm 



n 

1=1 



U\2a 



A'-k 



1/13 



X 



We note that ergodicity of A implies Ai > Xgeo > 1 (see (|30|) . (fT9|) and Proposition which ensures 
non-triviality of this lower bound. 

Now in order to finish the proof is suffices to combine the estimates (|38|) , (|48|) and (|53)) , and note 
that 



A 



2an 
geo 



^2a 



HA) 

d 



n 2ah(A)n 



which yields (|^ . 



3.3. Proofs of Theorem 121 ii) , iii) and Theorem |3j In this section we apply Theorem |3] to 
prove main theorems of Section |2j 

In order to determine the dissipation time of T^^a one has to determine the asymptotics of 
when n goes to infinity. According to formulas H17|) and H18() this problem reduces to problem H21|) 
solved in previous sections. 
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Thus in view of Theorem ^ there exist constants Ci and C2 such that for any 5,5' > and 
sufficiently large n 

Using formula (|17() 



e 

Now when n = n^j^^, we have 
and 

^—^ ( ln(l/e) - lnC2 I < Udiss < ^— 

(1 + 5')2ah{A) V / " " (1 - 6)2ah{A) 

which proves part ii) of Theorem [2 i-e. the logarithmic growth of dissipation time as a function of 
1/e. 

Moreover, using the definition of dissipation rate constant 

J-, IT'diss 

^diss — iim ■ , ^ 
e^o ln(l/e) 

we obtain 

1 



< Rdiss < 



(1 + 5')2ah{A) ~ ~ (1 - 5)2ah{A) 
Finally letting 5 — > and 5' ^ we arrive at the following results: 
• The general case - Theorem |21 iii) 

' <R.,.< ' 



2ah{F) 2ah{F) 
• The diagonalizable case - Theorem |31 

1 



diss 



2ah{F) 

This completes the proof. B 

4. Degenerate Noise 

In this section we compute the dissipation time for non-strictly contracting generalizations of 
a-stable transition operators. Instead of considering standard a-stable kernels of the form one 
can allow for some degree of degeneracy of noise in chosen directions by introducing the following 
family of noise kernels 

(55) 9e,a,B{^) ■■= Yl e~''''''''"ek(x), 

Where B denotes any d x d matrix with det B = 0. 

We denote by Ge^a^B the noise operator associated with ge,a,B- The degeneracy of B immediately 
implies that ||Ge^Q,^s|| = 1 and hence the general considerations of sections 1 and 2 do not apply 
here. The answer to the question whether or not the dissipation time is finite depends on the choice 
of matrix B. 

For simplicity we concentrate on the case when B is diagonalizable. 

We call the eigenvector of B nondegenerate if it corresponds to nonzero eigenvalue. 
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Theorem 5. Let F be any toral automorphism and T^^a,B = Ge,a,BUF- Assume that B is diago- 
nalizahle. Then 

i) If all nondegenerate eigenvectors of B* lie in one proper invariant subspace of F then dissi- 
pation does not take place i.e. ridiss = oo. 

ii) Otherwise the following statements hold. 

a) T^^a has simple dissipation time iff F is not ergodic. 

b) Ti^^a has logarithmic dissipation time iff F is ergodic. 

c) If T^^a,B has logarithmic dissipation time then the dissipation rate constant satisfies the fol- 
lowing bounds 

1 1 

^ Rdiss ^ 



2ah{F) ~ ~ 2ah{F) 
with some constant h[F) < h{F). The equality is achieved for all diagonalizable automorphisms F. 
Proof. 

We continue to use the convention A = F^ . The general formula derived previously for IITJJqII 
(see H17|) '). will now take the form 

(56) \\T^,^,b\\ = sup e-^S"-l^^''^l'"=e o^.«.^-^' 

Thus we need to estimate 

n 

inf VlSA'kP". 

To this end we denote by /ij {j = l,...,d) the eigenvalues of B and we construct a basis (of C^) 
{vj}j=i,...,d composed of normalized eigenvectors corresponding to eigenvalues /ij. We denote by 
Pjj=i d eigen-projections on Vj, and by Lj the set of corresponding linear forms, given 

by the eigenvectors Uj of B\ which are of course co-orthogonal to Vj, i.e. (uj, Vj) = for i ^ j. 
We have 

d d d 

X = ^ PjX = ^(Ljx)vj = ^(x, Uj)vj , Vx G W^. 

j=i j=i j=i 

In subsequent computations the symbols Ci, C2 denote some absolute constants values of which 
are subject to change during calculations. 
We consider two cases. 

i) All nondegenerate eigenvectors of B^ lie in one proper subspace of F. We have the following 
estimates 

d d d d 

\BA\\'^ > Ci^|P,-SA'k|2 = Ci^l^jf |P,-^'k|2 = C7i^|/xjf |(^'k,uj)|2 = Ci J^l^jf |(k,F'uj- 
j=i j=i j=i j=i 

and 

d d d d 

\BA'k\^ < C2Y,\PjBA'k\' = C2 J^l^.lVi^'kp = C2^|M,f K^'k,u,)|2 = C2 J] |^,pKk,F'u, 

j=i i=i i=i i=i 

Since at least one of ^ij is zero and all nondegenerate vectors Uj lie in a proper invariant subspace 
of F, one easily sees that for each fixed n 

n n d 

inf y|S^^k|2°= inf y y |u,f"|(k,F'u,-)|2" =0. 
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ii) In this case we have the following upper bound 



n n 



(57) inf y l^yl'kp" < ||5||2" inf V I A'kP" = ll^f" min V lA'kP". 

In order to provide an appropriate lower bound we note that the set of vectors {F^Uj}, where 
1 < h < d and j runs through the indices of all nondegenerate eigenvectors of B, spans the whole 
space (otherwise all nondegenerate uj would lie in one proper invariant subspace of F). We denote 
by {F'^'Uj.} (1 < i < d) a basis extracted from the above set. We can define now a new norm | • |u 
on R'^ by 

d 

|x|2 =^|(x,F'^-u,,)p 
1=1 

and compute 

dn n—1 d n—1 d d 

1=1 1=0 h=l 1=0 h=l 3=1 

n—1 d d n—1 d 
= Ci j;|/i,f"|P,-A^'+'^k|2- > j;|(^'^^+'^»k,U,-J|2" 

1=0 h=l j=l 1=0 i=l 
n—1 d n—1 

= Ci^^|(.4'^'k,F''»Uj-JP° = Ci^|A^'k|2" 

(=0 i=l 1=0 

Using the equivalence between norms | • | and | • |u and combining (|57)) with the above estimate we 
get 

n—1 dn dn 

Ci min Y\AV"< inf V Ifiyl'kP" < llsf " min V M'kP". 

This together with the obvious fact that h{A'^) = dh{A) and the general estimate (|22j) reduces the 
proof back to the nondegenerate case considered in the previous section. ■ 



5. Time of decay of ^-coarse-grained states 

The uncertainties in the initial preparation and the final measurement of the noiseless system 
give rise to non-cumulative random perturbations to the system. Alternatively, one can coarse-grain 
the initial and final states of the noiseless system by convoluting with the e- noise kernel. That is, 
instead of the original operator T^ (y^ we consider the operator Tg^Q, defined as 

(58) Tj^Q, = Ge^aUpGe^a = Gs^aUp'^Gg^a- 

and compute the number of iterations required to have the L^-norm of the final state being 
times that of the initial state. We will show that for ergodic toral automorphisms the required 
number of iterations is essentially the same asymptotically as the dissipation time computed in the 
previous sections. 

One can represent the action of Up or more generally Up in the Fourier series 

Upek = ^ ^^k^k'^k', 
o^k'ezd 
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where U]^\^i coincides with nk,k' defined previously (cf. <\l'S\i ) and 

"S' = Yl '*^k,ki'«ki,k2-"Wk,_i,k' 

0^ki,...,k„_igZ'' 

which satisfies 

(59) Yl l4"i'l' = l' Vn, k. 

o^k'eZ'* 

Then 

0^k„GZd 

Now we define 

SniK) = {ko G Z^{0} : ^xgk^ / 0}. 

Similar computations to these performed in Section [21 give the following general upper bound for 

1 1 1 1 



(60) iiT;vf< Yl E E 1^ 



W 1 2 

07^k„GZ<*koGS„{k„) koeS„(k„) 



For a toral automorphism one easily sees that 

(61) <;k„ = e-(l''«l^"+l^'''^°'^")5k„,A«ko 

and hence 



The arithmetic minimization problem (|18() corresponding to the dissipation time of T^^^ now be- 
comes 



(62) min (|k|2° + |A"k|2") . 

The key observation is that, by the same arguments as before, similar estimates to these given in 
(EH) hold 



(63) Ci I Y < |k|2° + |^"k|2° < C2 Af "|Ljk| 



2a 



The remaining computations are the same verbatim so the dissipation time of T^ q^ and T^^q^ are 
equal asymptotically. 

6. Time Scales in Kinematic Dynamo 

In this section we briefly discuss the connection between the dissipation time and some charac- 
teristic time scales associated with kinematic dynamo, which concerns the generation of electro- 
magnetic fields by mechanical motion. For a general setup and discussion we refer the reader to 
and ^S] and references therein. Here we restrict ourselves only to necessary definitions. 

Let B e Ll{T'^,W^) denote periodic, zero mean and divergence free magnetic field and let F be 
the time-1 map associated with the fiuid velocity. We define the push- forward map 

F,B(x) = dF(F-i(x))B(F-i(x)). 
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The noisy push-forward map P^^a on i^oC^'^'^'^) then given by 

(64) P,,«B := Ge,aF.-B, 

where the convohition (the action of Ge^a) is apphed component-wise. 

It is said that the kinematic dynamo action (positive dynamo effect) occurs if the dynamo growth 
rate is positive i.e. 

Rdyn = hm - In II P;"^ II > 0. 

n— >oo n 

Moreover if 

hm hm -ln||P" || > 0, 

then the dynamo action is said to be fast; otherwise it is slow. The anti-dynamo action takes place 
if 

lim -ln||P" II < 0. 

n—>oo n ' 

Now we introduce the threshold time scale as 

nth = max{n : ||P"q,II > ^ such that ||P"a^ll ll-Pe"^^ll ^ 

The threshold time n-this) is of order 0(1) as e — >■ for all fast kinematic dynamo systems. In 
the case of anti-dynamo action, rithis) captures the longest time scale on which the generation of 
the magnetic field still takes place. Finally nthi^) is not defined for systems with do not exhibit 
any growth of magnetic field throughout the evolution. In the case of anti-dynamo we consider the 
time scale on which the generation of the magnetic field achieves its maximal value 

Hp = mm{n : ||P"^|| = sup ||P™al|}- 

m 

which is called the peak time of the anti-dynamo action. 

Our next theorem establishes the relation between Up, nth and n^iss for toral automorphisms. 
Thus dF = F and 

Pe,aB = g,,a * P(B O p-^). 

Theorem 6. Let F be any toral automorphism. Then 

i) If F is nonergodic and has positive entropy then for all < £ < Rdiss^^PF the fast dynamo 
action takes place with dynamo growth rate satisfying 

Rdyn = InpF - sRdlss ^ InpF > 0, 

where pp denotes the spectral radius of F. The threshold time nth of order 0(1) and if F is 
diagonalizable then nth ~ [R-dyn^ ■ 

a) If F is nonergodic and has zero entropy then anti-dynamo action occurs and for nondiago- 
nalizable F, 

nth 



ln(nt/i) 

~ ndiss 

£ 

Moreover there exists a constant < 7 < d such that ||Pe"S|| ~ (1/e)'*'. If F is diagonalizable then 
||P"^|| is strictly decreasing (in n) and, hence, np = and nth is not defined. 
Hi) If F is ergodic then anti-dynamo action occurs and 



Hp ~ ndiss- 
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In particular if F is diagonalizable then 

np ~ nth - Rdiss In(ntfe) 
~ ^dissln(l/£) 

= ln(l/£) 

2ah{F) 

and 

(65) llPe^^ll - 

We see that even in the case of anti-dynamo action the magnetic field can still grow to relatively 
large magnitude when the noise is small (power-law in 1/e) . 

Proof. Representing the initial magnetic field B = (bi, ...,bd) in Fourier basis 

B= ^ B(k)ek 

one obtains 

P,,«B= ^B(k)e-I^'^l'"e^k, 
where we set A = {F~^)^. After n iterations 

Pe,a^= E i^"B(k)e-^Sr=il^'k|-g^„j^_ 

Thus 

II £,a II - ^ I V /I -QT^keZd ^ ' ^ ^' 

= g-2£min,_,^ggdEr=il^'k|2«|^nB|2 ^ g-2£ Er=i l^'k„|2" |^„g|2 < g-2£Er=i ||^n||2|B|2^ 

where k„ denotes a solution of the minimization problem 

n 

min VlA'kP". 
The above calculation provides the following upper bound 

IIP" II < p-£Er=il^'k„|2-||pn|| 
IP e,all — '-■ IK II- 

On the other hand let v„ denote a unit vector satisfying \\F^\\ = l^^vl. One immediately sees that 
the above upper bound for ||P"q,|| is achieved for magnetic field of the form B = v„ek„. Thus 

(66) IIP-JI =e-^Sr.ilA'k„r||^n||_ 

Now we consider the cases mentioned in the statement of the theorem 
i) Nonergodic, nonzero entropy case. 
For any nonergodic map we have 

n 

1=1 

This implies the following asymptotics 
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where ci,C2 > are constants (both equal iff F is diagonahzable). Thus for e < Rdiss^^PF we 
have 

Rdyn = In pF - eRdiss ^ In PF > 0. 
The threshold time is clearly of order 0(1) and in diagonahzable case can be written as 

1 e-^O 1 

In pp-sRj^^ InpF 

ii) Nonergodic, zero entropy case. 

In this case Inpp = 0. Thus if F is nondiagonalizable then (|67j) reads 
IIP" II PS e"^-^*^-^" I IF" 1 1 fsi e"^-^d»L"+'=il°"+c2 

II £,a II IK II ) 

with < ci < d. This immediately yields 

ci nth 1 



Rd: 



■ISS 1 



e ' \n{nth) e ' 
And moreover ~ (l/e)*^^. 

If F is diagonahzable then ||-F"|| = 1 and in this case ||-P"q,|| ~ e~^^diss'^ which implies Up = 0. 
iii) If F is diagonahzable, then form H22|l we know that for any < 6 < 1 and sufficiently large n 

n 

(68) yr_ {l-S)2ah(A)n ^ ^-^ y ^^l^^2a ^ ^{l+5)2ah{A)n ^ y^ 

Thus for large n we have 

maxe"^^+*p^ < max||P" || < max e"^'*'-VF- 

n n ' n 

We obtain the following constraints for Up 



lnA_|_5 \lnA_i_5y lnA_5 \e J ~ ~ InX^s \^nX_sJ lnA_5 Ve 
This gives 

1 , n„ 1 
< hm — < 



InA+a £-^oln(l/e) InA-^' 
Now since \±s — > for the above estimation yields the following asymptotics 

ln(l/e) PS Udissi nth - Rdiss ^n{nth) ~ Udiss- 



2ah{F) 

Similar asymptotic estimates (except for the constant) hold for nondiagonalizable F. ■ 

Appendices 

Appendix A. Affine transformations 

In this appendix we present a slight generalization of the results obtained in the paper. We 
consider here general affine transformations of the torus. The term affine transformations will 
be used here to refer to homeomorphisms of the torus with zero periodic but not necessary zero 
constant part (cf. Section I2.3|) i.e. transformations of the form F = F + c, where F is a toral 
automorphism and c is a constant shift vector. 

We begin with a short discussion of the ergodicity of affine transforms. 

The relation between ergodicity of a given affine transform F and associated with it toral auto- 
morphism F is summarized in the following proposition (for the proof we refer to appendix B) 
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Proposition 8. Let F be any toral automorphism. Then 

i) If F is ergodic then F is also ergodic. 

a) If F is not ergodic then F is ergodic iff 1 is the only root of unity in the spectrum of F and 
c • k Z'^ for any integer eigenvector k of F^ . 

Proof, i) Assume F is ergodic and for some c, F = F + c is not ergodic. Then there exists 
non-constant / G Lq(T'^) satisfying f = f o F ov in the Fourier representation 

(69) Yl = Yl e2"^^"'^"/(^~'k)ek 

where A = FK Comparing the absolute values of the coefficients we get 

(70) |/(k)| = |/(^-"k)| 

for any integer n and any k. However, ergodicity of F implies that A~"k ^ k for all k 7^ 0, which 
contradicts our assumption that / S Lq{T'^). 

ii) We will use the following fact, which can be proved by simple application of rational canonical 
decomposition. For any A E SL{d,'L) the following conditions are equivalent 

a) A possesses in its spectrum a root of unity not equal to one. 

b) There exists nonzero k £ Z'^ and a positive integer n such that k + Ak + ... + yl'^^^k = 0. 
Now assume that 1 is the only root of unity in spectrum of F (and hence of A) and c • k Z'^ 

for any integer eigenvector k of A, and that both F and F are not ergodic. The latter assumption 
implies the existence of a non-constant / E -Lq(T'^) satisfying equations (|69|) and (|70j) . Relation 
(|70|) clearly imphes that if /(k) / then >l"k = k for some n. Moreover, since 1 is the only root of 
unity in spectrum of j4, we have, in view of b) that Ak = k. Thus the only possible non-constant 
invariant functions of F are single Fourier modes corresponding to integer eigenvectors of A. But 
if such a Fourier mode is invariant under F then directly form (|69|) one concludes that e^'^*'*"'^ = 1 
or equivalently k - c € Z*^, for some integer eigenvector of A. To prove the converse we assume that 
F is not ergodic and consider two cases: 

Case 1. A possesses in its spectrum a root of unity not equal to one. In this case according to 
condition b) there exists nonzero k G Z'^ and a positive integer n such that k + ylk+... + 74"'^"'^k = 0, 
which implies in particular that ^"k = k and Ak 7^ k. Now we define the function 

J=ek + e BAk + .-. + e \^i-o ) e^n-ik 

which clearly satisfies the condition / = f o F. This proves that F is not ergodic. 

Case 2. There exists integer eigenvector of A such that k • c G Z'^. Then clearly for such k, 
/ = ek is F-invariant and hence again F is not ergodic. ■ 

We recall that c = (ci,..,Crf) generates ergodic shift on the torus iff l,ci,..,Cd are linearly 
independent over rationals. Thus as a direct consequence of the above proposition we get 

Corollary 5. If F is not ergodic and 1 is the only root of unity in the spectrum of F then F is 
ergodic for all vectors c generating ergodic shifts on the torus. 

Now we are in a position to state and prove the generalization of Theorem |21 from Section [2.31 to 
the case of affine transforms (the corresponding generalizations of Theorem and Corollary are 
straightforward) 

Theorem 7. Let F be any affine transformation on the torus T*^, F associated with F toral auto- 
morphism and TfT^Q, = Gs^aUp. Then 

i) Tfr^a has simple dissipation time iff F is not ergodic. 

'i-'i) Te^a has logarithmic dissipation time iff F is ergodic. 
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Hi) IfT^^a has logarithmic dissipation time then the dissipation rate constant satisfies the follow- 
ing constraint 



2ah{F) 2ah{F) ' 

where h{F) < h{F) is certain positive constant. 

Remark 1. The dissipation time of an affine transformation F is determined by ergodic properties 
of its linear part F and hence not by ergodic properties of F itself. In particular all ergodic affine 
transformations associated with nonergodic toral automorphisms ( cf. Proposition 0) have simple 
dissipation time. 

Proof of Theorem [7| Specializing the general calculations of dissipation time presented in 
Section 12.21 to the case of affine transformations F = F + c, with nonzero c, one easily finds the 
following counterparts of formulas ((TB|) and H17() 

_27rik-cr 

t^„(ko,k„j — e ) e ^'-i ' ' OA"ko,k„- 

Now, in order to determine the dissipation time of T^^a = Ge^aUp one has to determine the asymp- 
totics of llT^all as n goes to infinity. According to the above formulas and formulas ()15|) and (jlZj) 
from Section the value of HTJ^all does not depend on c, which reduces the proof to the case 
c = considered in the main body of the paper. ■ 

Appendix B. Proofs of some elementary facts 
Proof of Proposition [2] 

The proof will be based on the Riesz convexity theorem (see HTj? PP- 93-100) which states that 
for any operator T defined on L*'(T'^), 1 < p < cxd, In ||T||p is a convex function of On the space 
LP(T'^) we consider the operator T := T^ ^ — (") and we have the relation T"/ = T^^^if ~ {f))yf ^ 
LP{T'^),n > 1 because T^^a is conservative. Now since ||/ — {f)\\p < 2||/||p, it follows that 

(71) llT'^llp < 2||T;^„||p,o < 2 

(72) wT^jpfi < \n 

for l<p<oo,n>l. The Riesz convexity theorem implies that if p < g < oo 

(73) In ||f "II, < ^ In ||r"||p In ||r"||oo 



while if 1 < g < p 

(74) In ||f "II, < (i^) In ||f "||, + (l - In ||f "||,. 

^From H73|) -(|74 |) we have the interpolation relations 

(75) ||f"||, < ||r"||P/5||f"||^P/5, p<q<oo 

(76) ||f"||, < ||r"||(i-^"')/(i-P"')||r"||}-('-^")/('~P"'\ l<q<p 
which, along with ((7T|) - ((7^ . imply 

\Ka\ko < nT^Xl p<q<^ 

< 2\\T-X,'^-'^/^'-^''\ l<q<p 
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O^kiGZd 



This proves that the order of divergence of Udissip) are the same for 1 < p < oo. Estimates 
((75|) - (|7H|) also show that the order of divergence of n^j^sCl) and n-iHssioo) is at least as high as 
n-dissip), 1 < p < oo. ■ 

Proof of Lemma ^ 

Using the notation introduced in Section 12.21 one has 

^^ko,kiWki,k2---'«k„_i,k„e~^^"=il'''l'"ek„ = Yl ^n(ko,k„)ek„. 

07^ki,...,k„GZ'« 0^k„GZ'* 

We note that for any n and k„ € Z*^, the sequence W„(ko,k„) (indexed by ko € Z*^) belongs to 
/^(Z"^). Indeed, using the Cauchy-Schwarz inequality and identity H13() one gets for n = 2, 

2 

Y |^^2(ko,k2)P= Y 

O^ko GZ*^ OT^ki GZ'* OT^ki GZ"* 

< g-e|ki|2" ^ ^_e|kip-g-2e|k2|2" ^ ^g-2£|k2p-^ 

07^kiGZ<^ O^kiGZ'* 

where K denotes a constant. Similar estimates hold for n > 2. 
Now applying the Cauchy-Schwarz inequality in (|14|) we get 

(77) \\T^,jf < E E l/>o)l' E K(ko,k„)|^ 

OT^k^GZ-^koG^nCkn) koG5„(k„) 

■ 

Proof of part i) of Theorem ^ 

In view of theorem ^ it suffices to construct an eigenfunction of Up which belongs to Lq{T'^) n 
H^a^jdy Directly from Proposition El one concludes that F, and hence also A, possesses a root 
of unity in its spectrum. This means that j4™'ko = ko, for some m and certain nonzero vector ko, 
which can be chosen to be an integer. Now we define 

/ = eko + e^ko + ••• + e^^-iko 
Obviously / G Lq(T'^) n //^"(T'^), for any a. To complete the proof it suffices to notice that 

Upf = e^ko + e^ako + ••• + ^A^ko = ^ko + e^ko + ••• + e^^-iko ^ f- 



Proof of Proposition |3J 

For the purposes of the proof we use the following abbreviation 

• PRSiW^) - proper rational subspace of M'^. 

• PIS{A,W^) - proper ^-invariant subspace of M"^. 

• PRIS{A,W^) - proper rational ^-invariant subspace of M'^. 

a) ^ b): . Suppose there exists PRIS{A,M.'^) Si. Let Ai be a matrix representing invari- 
ant rational block associated with Si. Then Ai is rational matrix and its characteristic 
polynomial Pi belongs to Q[x]. Let P denote the characteristic polynomial of A. Then 
P = P1P2 and since both P, Pi G Q[x] then also P2 G Q[a^], which means P and hence A is 
not irreducible. 
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b) ^ c): . Assume there exists PRS{W^) S contained in PIS{A,M.'^) V. Take any rational 
vector q G 5 and let do = dimV then the set {q,^q, .., A'^O'^q} spans PRIS{A,W'-). 

c) ^ d): . Assume that for given q and an arithmetic sequence ni,...,n[i, the set S = 
{^"iq, A^^q, ^"''q} does not form a basis. Since for some fixed integer r, ni = ni + (l — 
l)r, we have A^'q = {A'^y~^A'^^q = (A^)'~^q, where q = ^"^q. Now consider the biggest 
subset So = {q, A^q, (A'^)^q, (^'')'^°~^q} such that do < d and 5*0 is linearly independent. 
Obviously So spans a PRIS{A'') which is also a PRIS{A). 

d) ^ a): . Suppose that characteristic polynomial P of ^ is not irreducible in Q[x]. Then 
P = P1P2, with Pi,P2 G Q.[x]. From the Cayley-Hamilton theorem we get that = 
P{A) = Pi{A)P2{A). Hence for any nonzero rational vector q, either 1) P2{A)c[ = or 2) 
q := P2{A)c[ / and Pi(yl)q = 0. Since max{deg{Pi, P2)} < d, there exists a nonzero 
rational vector q (namely q or q) such that the set of iterations {q, ^q, A^q, ^'^^-'^q} 
does not form a basis of M.'^. 

e) ^ f): . Assume there exist nonzero q G Q'^ orthogonal to certain PIS{A,M.'^) V. Then for 
any n and any / e V, ((At)"q,/) = (q,^"/) = and hence S = {q, ^l+q, (At)2q, (^t)'i-iq}^ 
cannot form a basis, which in view of equivalence a) <^ d) implies reducibility of . 

f) g): . Suppose there exists PIS{A,'R.'^) V contained in certain PRS{W^) S. Since S is 
rational, S-^ is also rational. Consider any rational vector q G S^, then (q, /) = for any 

g) =^ b): . If there exists PRIS{W^), then this subspace is yl-invariant and contained in 
PRS{R'^) i.e in itself. 

Now since b) is equivalent to a) it is enough to establish the equivalence between a) and e) to 
complete the proof. But the latter equivalence is obvious in view of the fact that A and A^ have 
the same characteristic polynomial. ■ 

Proof of Proposition \E\ 

Suppose ^4 is a toral automorphism of zero entropy. The latter property is equivalent to the 
fact that all the eigenvalues of A are of modulus 1. Let Pa be a characteristic polynomial of A. 
Consider any irreducible over Z factor P of polynomial Pa and construct a toral automorphism B 
such that its characteristic polynomial is equal to P . Obviously all the eigenvalues of B are also 
the eigenvalues of A, and each eigenvalue of A can be found among eigenvalues of some matrix B 
of this type. Irreducibility of P implies irreducibility and hence diagonalizability of B. 

Thus for any nonzero vector k S Z'^ and any positive integer n the following estimate holds 
|i?"k| < |k|, which implies the existence (for each k) of some integer r such that B^"k = k. 

The latter shows that all the eigenvalues of B (and hence also of A) are roots of unity. ■ 

Proof of Proposition El 

We first show that irreducible polynomials P G Q[x] do not have repeated roots. Indeed suppose 
A is a root of P of multiplicity greater that 1, then A is also a root of a derivative polynomial 
P' S Q[x]. Since the minimal polynomial of A must divide both P and P' and deg{P') < deg{P) 
one immediately concludes that P is not irreducible. Now, suppose A S GL{d, Q) is completely 
decomposable over Q and let (|25|) be its block diagonal decomposition into irreducible blocks. 
Each Pa^, as a characteristic polynomial of Aj, is irreducible over Q and hence does not possesses 
repeated roots, which implies diagonalizability of each Aj and hence of A. On the other hand if A 
is diagonalizable then its minimal polynomial does not possesses repeated roots, which implies that 
all characteristic polynomials associated with elementary divisors are (first powers of) irreducible 
polynomials. This implies irreducibility of each block in representation H25() . ■ 

Proof of Proposition [Tj 
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Let Pa be the characteristic polynomial of an irreducible matrix A E GL(d,Q). Since Pa is an 
irreducible element of Q[x] it does not possesses repeated roots (see the proof of Proposition IHl). 
■ 
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